Extracting the Mott gap from energy measurements in trapped atomic gases 
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We show that the measure of the so-called release energy, which is an experimentally accessible 
quantity, makes it possible to assess the value of the Mott gap in the presence of the confinement 
potential that is unavoidable in the actual experimental setup. Indeed, the curve of the release 
energy as a function of the total number of particles shows kinks that are directly related to the 
existence of excitation gaps. Calculations are presented within the Gutzwiller approach, but the 
final results go beyond this simple approximation and represent a genuine feature of the real system. 
In the case of harmonic confinement, the Mott gaps may be renormalized with respect to the uniform 
Ccise. On the other hand, in the case of the recently proposed off-diagonal confinement, our results 
show good agreement with the homogeneous case. 
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I. INTRODUCTION 

It has been established that ultra-cold Bose and Fermi 
gases trapped in optical lattices provide experimental re- 
alizations of long-standing lattice models widely consid- 
ered in condensed matter physics and statistical mechan- 
ics, such as the Bose and Fermi Hubbard models [l|, Hi- 
One of the most spectacular achievements was the ob- 
servation of the superfluid to Mott-insulator transition 
as the optical lattice depth is varied Q. This transition 
has been explained in light of the Bose-Hubbard model 
with on-site repulsive interactions and hopping between 
nearest-neighboring sites. In the superfluid phase, each 
atom is spread out over the entire lattice, with long-range 
phase coherence, whereas in the insulating phase, exact 
numbers of atoms are localized at individual lattice sites, 
with no phase coherence across the lattice; this phase is 
characterized by a gap in the excitation spectrum and a 
vanishing compressibility. Similarly, recent experiments 
have provided evidence on the formation of a Mott insu- 
lator of fermionic atoms in an optical lattice [j] . Several 
challenges arise when quantitatively comparing experi- 
mental data and theoretical results obtained from the 
analysis of infinite systems. One of the most important 
problems is due to inevitable spatial inhomogeneities in- 
duced by the optical trap, which is necessary to confine 
particles la-Q- Further complications are connected to 
determine temperature effects present in experiments |8| , 
and also to the limited available tools for the experimen- 
tal characterization of the phases. 

Here, we would like to focus our attention on the pos- 
sibility to make quantitative estimations of the Mott 
gap Eg from relatively simple quantities that can be 
experimentally addressed. In presence of the optical 
trap, it is usually not possible to have a Mott insulat- 
ing phase throughout the whole lattice and compress- 
ible regions intrude the system. This fact has been 
widely discussed both experimentally and theoretically 
and shows up through the typical "wedding-cake" pro- 
file of density 0, [J0, i, ^. Therefore, the system 



always possesses regions that are locally compressible 
and a precise determination of the Mott gap is subtler 
than in the homogeneous case. Some preliminary at- 
tempts to measure the excitation spectrum of interact- 
ing bosons have been performed by using Bragg spec- 
troscopy [ll|. Other approaches to characterize the ap- 
pearance of Mott-insulating regions from experimen tally 
accessible quantities have been also proposed [l2|, [l3j . 
Here, we would like to propose an alternative approach 
that is based upon energy measurements only and could 
give important insights about the appearance of Mott- 
insulating regions, as well as the actual value of the gap. 

For an infinite and homogeneous system, the excita- 
tion gap Eg can be calculated from the knowledge of 
the total energy for different particle numbers, namely 
Eg ^ n^ — ji^ , where fj,^ = ±(i57v±i — E'at), En being 
the ground-state energy with N particles. The Mott gap 
Eg is finite whenever /i+ 7^ fi^ and, therefore, introduces 
a discontinuity in the first derivative of the energy with 
respect to the density. 

In this paper, we will show that the existence of in- 
compressible Mott regions and the values of the corre- 
sponding gaps can be obtained from so-called release en- 
ergy E'^'^'', which may be measured experimentally (l4| . 
Indeed, E'^'^^ is obtained by integrating the momentum 
distribution of the atoms after having switched off the 
confinement and let the atoms expand freely. The release 
energy is given by the sum of the kinetic and interaction 
energies just before switching off the trap |14| 
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We notice that it would be much more difficult to extract 
the total energy, E*"* ^ E''"' + E"''* + EP°^ that also 
includes the potential term due to the trap, since this 
would require the knowledge of the density profile in the 
presence of the trap, which is hard to reconstruct. 

Specifically, we address two types of confinement i) the 
usual harmonic confinement and ii) a recently proposed 
off-diagonal confinement [l5|. In such a confinement, the 
strength of the hopping parameter is varied across the 



lattice, being maximum at the center of the lattice and 
vanishing at its edges, which naturally induces a trapping 
of the particles. 

We consider bosons loaded in one- and two- 
dimensional optical lattices and use an insightful vari- 
ational approach based upon the Gutzwiller wave func- 
tion [la . llTJ I . Similar results must hold also in fermionic 
systems. Moreover, although approximated, this varia- 
tional method is expected to correctly capture the behav- 
ior of the exact ground state. We find that the presence 
of Mott regions in the system is signaled by discontinu- 
ities in the derivative of the release-energy curve with 
respect to the total number of bosons, reminiscent of the 
presence of a gap in infinite homogeneous system. In the 
case of harmonic confinement, the measured gap may be 
substantially smaller than the one of the uniform system, 
whereas a much closer agreement is achieved by consid- 
ering the off-diagonal confinement. 

The emergence of strong signatures due to the forma- 
tion of Mott domains in the energy measurements has 
been explored in Ref. |13| . where it has been shown that 
the appearance of Mott domains is accompanied by min- 
ima in the release-energy curve as function of the on-site 
interaction parameter. Here, we make a further step in 
this direction and provide a direct and quantitative con- 
nection between energy calculations and Mott gaps, in 
presence of an external trap. 

The paper is organized as follows: in section |lll we in- 
troduce the Bose-Hubbard model and describe the varia- 
tional method, in section Hill we present our results, and 
finally, in section ITVl we draw our conclusions. 



II. MODEL AND METHOD 

Our starting point is the Bose-Hubbard model which 
describes interacting bosons on a lattice [6[: 
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where (...) indicates nearest-neighbor sites, bj (bi) cre- 
ates (destroys) a boson on site i, and hi is the local den- 
sity operator. U is the on-site interaction, tij is the 
hopping amplitude, and e^ is a local energy offset due to 
an external trapping potential. 

To study this model and describe its ground-state 
properties we use a mean-field approximation based on 
the Gutzwiller ansatz [16|, Il3l ■ This simple approach is 
able to capture important features of the true ground 
state and provides a qualitatively correct description of 
quantities such as the local density or the total energy, 
even in presence of spatial inhomogeneities |6Jil10i 16. 17|. 
More involved wave functions with a long-range Jastrow 
factor [l8| , or numerically exact calculations [l9[ may be 
considered, but a much larger computational effort would 
be required. Instead, here we are just interested in show- 
ing that Mott gaps can be extracted from the behavior 




FIG. 1: (Color on-line) Energy per site versus density in the 
case of tlie ID homogeneous model for different values of U/t. 
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FIG. 2: (Color on-line) Upper panel: Release energy. Middle 
Panel: local density Ui across the lattice sites. Lower Panel: 
local density fluctuations Ai across the lattice sites. A value 
of Vo/t = 0.01 is considered. AH quantities are shown as a 
function of the total number of bosons A'' and the interaction 
is set to U/t — 4. Calculations are shown for the ID model. 



of the release energy as a function of the total number 
of particles, and we do not expect qualitative differences 
when considering more accurate approaches. 

Within the Gutzwiller ansatz the ground-state wave 
function is approximated as 



\^o) = u j:.a 
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where Ito)^ is the Fock state with m particles at site 
i and /^ are variational parameters which have to be 
determined by miniinizing the expectation value of the 
Gutzwiller ansatz on the Hamiltonian in Eq. ^. The 
sum in Eq. ([3|) runs from states with zero particles up to 
infinity; however, from a numerical point of view, we have 
to consider a cutoff and take only states up to a maximum 
number of particles (per site) M^^^ ^ (^gI^^iI^g) such 



that the contribution of those states with higher density 
are neghgible and observables are converged to a certain 
desired precision. 

Equivalently, the Gutzwiller wave function can be in- 
troduced as the ground state of the following mean-field 
Hamihonian [20,1^: 



H 



mf 



\ E ^^.^ (^I*^ + **^^ - ***^) + ^■''■ 



(»-j> 



u 



'^ni{ni-l) + '^{ei- ^l)n^, (4) 



where ^^ is the mean-field potential which is self- 
consistently defined as *i = {^c\bi\^G)] it can be shown 
that /*^ is related to the ground-state eigenvector com- 
ponents of the converged solution of the local Hamilto- 
nian (0]) [2^. The parameter /i is the chemical potential 
that fixes the number of bosons. In the recent past, the 
Gutzwiller approach has been successfully used to study 
a wide range of bosonic systems, like for instance dis- 
ordered potentials [22[, boson-boson mixtures [23], and 
even time-dependent problems [2J-[26[. 



For the sake of simplicity and to simplify the presenta- 
tion of the results, we first consider the one-dimensional 
(ID) case, which is the limiting case where a collection 
of non-interacting tubes is created. An almost ID model 
may be easily generated by using different lasers in the 
three different spatial directions and has been experi- 
mentally considered [111, l23|- We mention that in ID, 
within the Gutzwiller approach, the compressible phase 
has inevitably a finite condensate fraction. However, this 
approach correctly reproduces the opening of the excita- 
tion gap at the Mott transition; therefore, the presence 
of a condensate does not imply relevant qualitative dif- 
ferences on the estimation of the gap with respect to an 
unbiased calculation. Furthermore, we also report some 
results for two dimensions (2D) for which we obtain simi- 
lar conclusions. We emphasize that all the results remain 
valid in higher spatial dimensions because the superfluid 
to Mott-insulator transition occurs in any dimension ac- 
companied by the opening of a gap in the spectrum. 

We consider a lattice with a harmonic potential of the 
form ei = Vorf, where r^ is the distance of site i from the 
center of the lattice. In this case, the hopping amplitude 
is kept constant for all lattice sites tij = t. In addition, 
we also analyze the case of off-diagonal confinement in 
which Ei ~ and t^.^+i = 4t x i x {L — i)/L'^, where L 
is the total number of sites. In both cases, we evaluate 
local quantities which allow us to determine whether a 
certain region across the lattice is in a compressible or 
incompressible state, as well as the release energy of the 
system. In particular, we will show, besides the release 
energy, the local density rii = ('I'G|ni|^G) and its fluctu- 
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FIG. 3: (Color on-line) The same as in Fig.[2]for U/t = 15. 
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FIG. 4: (Color on-line) Upper panels: first Mott gap obtained 
from the release energy as function of U/t for Vo/t = 0.01, 
for comparison the homogeneous case is also reported (left); 
release-energy curves for different U/t and Vo/t = 0.01. Lower 
Panel: first Mott gap as function of Vo/t for U/t = 15. Cal- 
culations are shown for the ID model. 



III. RESULTS 

Before showing the results for the confined system, 
which is relevant for experiments, let us briefly discuss 
the homogeneous case, with Cj = and tij = t. In 
this case, a superfluid-Mott transition takes place at in- 
teger fillings whenever the on-site interaction U is large 
enough. On the other hand, for any non-integer fillings 
the ground state is always superfluid and, therefore, com- 
pressible. In Fig. [U we report the energy curve as a 
function of the density n. Within the Gutzwiller ap- 
proximation the values of the critical interaction may be 
determined analytically, i.e., Uc/t = D[^Jn + \Jn + 1)^, 
where n is an integer [l^l. Whenever U < Uc, the en- 
ergy curve is smooth with a positive curvature, implying 
a finite compressibility and a vanishing gap. On the con- 
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FIG. 5: (Color on-line) Upper panels; 2D results on local 
density fluctuations Ai across the lattice before (left) and 
after (right) the kink in the release energy for U /t = 20. 
Lower panel: Release energy for different values of U /t and 
fixed Vo/t = 0.01. 



trary, for U > Uc, there is a discontinuity in the curve 
at integer fillings, (the behavior in the vicinity of 7i = 1 
is reported in Fig. [1]) , signaling the presence of the Mott 
gap. The latter one can be estimated by considering the 
change of the slope close to the discontinuity. 

Let us now switch on the harmonic confinenient by set- 
ting Vo/t = 0.01. In Figs. [2] and [3l we present the results 
for the release energy as function of the total number 
of bosons N, as well as the local quantities rii and A^ 
across the lattice sites. For the case with U/t = 4 (see 
Fig. [2]) there are no insulating phases, no matter what 
the number of bosons is. The density profile is smooth, 
with a broad maximum at the center of the trap. In this 
case, all regions of the lattice are (locally) compressible 
and, therefore, the ground state is gapless. This is not 
the case for U/t = 15 (see Fig. [3]), where insulating re- 
gions are expected. Indeed, what is found is the usual 
"wedding-cake" structure in the local density: the Mott 
regions with integer rii and vanishing A^ are surrounded 
by compressible regions which are locally gapless. We 
emphasize that the vanishing of A^ is consequence of the 
Gutzwiller approach; in a more accurate description of a 
Mott insulator, this quantity is indeed finite, though it is 
strongly reduced with respect to its value in the super- 
fluid regions. Here, the derivative of the release- energy 
curve clearly exhibits discontinuities that are reminiscent 
of the presence of a Mott gap. This discontinuity takes 
place whenever a new compressible region appears at the 
center of the trap, on top of the underlying Mott phase, 
see Fig.|3l In this way, we can define an energy gap for the 
confined system exactly as in the homogeneous system. 



namely i?-' = /x+, - M" „ where ^,f^, = ±{E^^^i, - E^^) 



(E"]^' being the release-energy with N particles). We 
would like to stress that similar results may be obtained 
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FIG. 6: (Color on-line) Results on total energy of the confined 
system in the 2D model as function of the number of particles 
for different values of U/t. 



can be considered in experiments. For the case shown in 
Fig. [3] (i.e., for Vo/t = 0.01 and U/t = 15), we obtain 
that the first Mott gap is E^'^'- /t ~ 7.2, to be compared 
with the Mott gap with n = 1 oi the homogeneous case 
that gives Eg/t ~ 11.8. The reduction of the measured 
gap comes from the fact that the release energy in the 
presence of the trap contains not only the information 
about the local creation of the new compressible region 
at the center of the trap, but also about all other sites 
of the lattice, which do not undergo the Mott transition. 
Therefore, the effect is spatially averaged over regions 
that are locally compressible and incompressible. Never- 
theless, the effect that originates from the central sites is 
visible and allows us to provide an estimate of the gap 
associated with such a transition. 

A summary of the results is reported in Fig. 21 where 
we show Eg'^'' (for the first Mott gap) as a function of U/t 
for Vo/t = 0.01 and as a function of Vo/t for U/t = 15. 
We find that, when the harmonic trap is increased, the 
Mott gap saturates to E'^J^^ /t ~ 10.6, which is close to the 
value of the homogeneous system. The initial depletion of 
the Mott gap as a function of Vb is due to the presence of 
(large) regions of compressible sites close to the borders of 
the system; by further increasing Vq, these regions shrink 
and the gap eventually tends to approach the value of the 
uniform case. 

We now briefiy consider the two-dimensional case, see 
Fig. [5] Exactly as in ID, in the appearance of a new com- 
pressible region at the center of the trap is accompanied 
by a discontinuity in the derivative of the release-energy 
curve. Similar conclusions to those in ID are obtained, 
confirming the appearance of kinks in the release-energy 
curve in 2D where the Gutzwiller mean-field approach is 
more reliable. 

We stress that in presence of a diagonal confinement 
is such that compressible and incompressible regions co- 
exist. Indeed, the true gap of the overall confined sys- 
tem will be zero, as the global compressibility is always 
finite [a, [12] ■ Indeed, the total energy, as opposed to 
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ties n = N/L, as in a homogeneous case. The advantage 
of this kind of confinement is that it is possible to have 
a Mott-insulating pha se throughout the whole lattice, as 
discussed in Ref. [l5|. By using the energy curve, we 
find that the value of the gap at n = 1 is E'^J'^ /t ~ 11.9, 
remarkably close to the corresponding value in the homo- 
geneous system {Eg/t ~ 11.8) but slightly larger. This 
is due to the fact that regions towards the border of 
the lattice are effectively in a deeper Mott phase due to 
their reduced kinetic energy. Away from integer densi- 
ties, the system shows coexistence of compressible and 
incompressible states as in the harmonic confinement, 
with compressible sites always at the center of the trap 
and Mott sites on the borders, again due to the reduced 
kinetic energy there. 



FIG. 7: (Color on-line) The same as in Fig. [3] but for an off- 
diagonal confinement. The number of sites is L = 200 and 
n = N/L is the density. 



IV. CONCLUSIONS 



the release energy, is a completely smooth curve with 
positive curvature (hence a vanishing real gap and finite 
compressibility). The total-energy curve as function of 
the total number of particles is shown in Fig. [6] for a 
2D system. This curve is completely smooth, even when 
Mott insulating domains are present in the trap [5|. 

Now, let us analyze the case of off-diagonal confinc- 
mentjTecently proposed by Rousseau and collabora- 
tors [13. In the experimental setup, the release- energy 
measurement corresponds to the actual total energy of 
the trapped system, since in this case there is no po- 
tential energy. In Fig. [3 we present the results for 
U/t = 15. We observe that the kinks in the energy curve 
occur exactly at values of the total number of bosons N 
which arc commensurate with L, i.e., for integer dcnsi- 



We have shown that relatively simple energy measure- 
ments may give important insights into the actual value 
of the Mott gap. In particular, it has been shown that 
the presence of the harmonic potential inay renormalize 
the value of the gap with respect to the uniform case. 
On the other hand, much closer results may be achieved 
when the off-diagonal confinement will be experimentally 
realized. Although the experimental determination of 
the Mott gap from release-energy curves may be diffi- 
cult, since small error bar on both energies and particle 
numbers are required, we hope that in the near future 
this method will be successfully applied. 

We acknowledge useful discussions with M. Fabrizio, 
L. FaUani, C. Fort, S. Giorgini, V. Rousseau, and A. 
Trombcttoni. 
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